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The model assumptions of the recently formulated framework of highly- 
anisotropic and strongly-dissipative hydrodynamics (ADHYDRO) are ana- 
lyzed. In particular, we study dependence of numerical results on different 
forms of the entropy source and compare our approach with other frame- 
works describing locally anisotropic fluids. We also discuss the effects of 
different forms of the initial conditions on the process of isotropization. 
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1. Introduction 

Recently, we have introduced the framework of highly-anisotropic and 
strongly-dissipative hydrodynamics [HE] (ADHYDRO). This framework is 
suitable for modeling of early stages of relativistic heavy-ion collisions, in the 
situations where very large pressure anisotropies are expected [31 HI [5] . In 
Ref. [1] our approach was used to describe boost-invariant, one-dimensional 
systems, whereas in Ref. \2\ we applied ADHYDRO in non-boost-invariant 
situations. 

ADHYDRO is a simple extension of the perfect-fluid hydrodynamics. 
For locally isotropic fluids it is reduced to perfect-fluid description, and 
for small pressure anisotropies and purely longitudinal expansion it agrees 
with the 2nd order Israel-Stewart viscous hydrodynamics (we note that 
longitudinal expansion dominates at the early stages of heavy- ion collisions). 
For highly-anisotropic systems, ADHYDRO may be treated as an effective 
model which describes strongly non-equilibrium dynamics. 

* Supported in part by the Polish Ministry of Science and Higher Education, grants 
N N202 288638 and N N202 263438. 
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The construction of ADHYDRO was motivated by the idea of develop- 
ing an approach that is suitable for uniform description of several stages 
of the evolution of matter — starting from the very early stages where the 
momentum anisotropics are extremely high, through the transient viscous 
evolution, and ending with the perfect-fluid stage. Since most of the real- 
istic hydrodynamic calculations indicate that the evolution of the system is 
very well described by the perfect-fluid hydrodynamics (or, equivalently, by 
viscous hydrodynamics with a very small ratio of the shear viscosity to en- 
tropy) with the starting proper time below Thyj = 1 f m 
we expect that the ADHYDRO description is reduced or becomes very close 
to perfect-fluid hydrodynamics at the proper time r w Thyd- On the other 
hand, at the very early stages of the collisions, r < Thyd, neither perfect- 
fluid hydrodynamics nor viscous hydrodynamics are formally applicable. In 
particular, viscous corrections become compatible with the leading terms, 
which may lead to negative pressure. Thus, the very early dynamics requires 
the use of the kinetic theory (which has its own well known limitations) or 
modeling. ADHYDRO is a proposal for such modeling. 

A work along similar lines have been recently done by Martinez and 
Strickland \12\ I13j . who derive analogous equations from the Boltzmann 
equation with the collision term treated in the relaxation time approxima- 
tion. Instead of the usual expansion around the isotropic one-particle dis- 
tribution, they consider an expansion around an anisotropic distribution. 
In this paper we compare some aspects of our approach with that presented 
in [UllIS]. 

The main model assumptions of ADHYDRO are the form of the energy- 
momentum tensor that allows for anisotropic pressure and the form of the 
entropy source. For small anisotropics, the form of the entropy source is 
restricted by the consistency with the 2nd order viscous hydrodynamics. On 
the other hand, for large anisotropics, various forms of the entropy source 
are conceivable. Of course it would be very useful to obtain some hints about 
the form of the entropy source from the underlying microscopic dynamics. 
Such connections are, however, difficult to find since the initial dynamics 
of the system produced in heavy-ion collisions is very much complicated. 
Therefore, in this paper we simply use different forms of the entropy source 
and study their effects on the evolution of the system. 

Besides the analysis of the entropy source term, we discuss the micro- 
scopic interpretation of our framework and analyze the effects of quantum 
statistics. We find that the inclusion of the Bose-Einstein statistics leads to 
minor quantitative modifications of the approach introduced in [T]. More- 
over, in this paper we present new results obtained for the initial pressure 



^ We use the natural units where h — c — 1. 
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anisotropy that depends on rapidity. We demonstrate that such initial con- 
ditions do not affect the isotropization process. 

2. Basic equations 

The ADHYDRO model is based on the following form of the energy- 
momentum tensor jTl| [T5] 

T^"' = {e + P±)U''U'' - P±g'"' - {P±- Pii)^^. (1) 

In Eq. ([1]) e, Pi, and Py are the energy density, transverse pressure, and 
longitudinal pressure, respectively. For Pi = P|| = P we recover the form 
valid for the perfect fluid. The four-vector describes the hydrodynamic 
flow 

U''=jil,V„Vy,V,), 7= (l-t;2)-V2^ (2) 

while defines the longitudinal direction (corresponding to the collision 
axis) 

^'^ = 7. (^.,0,0,1), 7. = (l-t'')"'/'- (3) 
The four-vectors and satisfy the following normalization conditions: 

U^ = l, V^ = -l, U-V = 0. (4) 

In the local-rest-frame (LRF) of the fluid element the four- vectors and 
have simple forms, 

^7^^ = (1,0,0,0), = (0,0,0,1), (5) 

and the energy- momentum tensor has the following diagonal structure 
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(6) 



In addition to the energy-momentum tensor ([T]) we introduce the entropy 
flux 

= aU>', (7) 

^ The Lorentz transformation leading to LRF is defined explicitly in Sect. 1111 Since 
this transformation should not change the longitudinal direction, it must be defined 
in a special way. 
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where a is the entropy density. In our approach e and a are functions of 
the two pressures, P± and Py. In particular, for massless partons we have 

e = 2P^ + Pll. (8) 

The space-time dynamics of the system is described by the fohowing equa- 
tions 

d^Tf'" = 0, (9) 
S^a'^ = (10) 

Equation ([9]) expresses the energy-momentum conservation, while Eq. (|10p 
describes the entropy production. The form of the entropy source S deter- 
mines the dynamics of the anisotropic system 0. This issue will be discussed 
in more detail below. 

If the functions (t(Pi,P||) and S(Pi,P||) are specified, Eqs. ([9]) and 
pop form a closed system of 5 equations for 5 unknown functions: three 
components of the fluid velocity v, the transverse pressure P±, and the 
longitudinal pressure The projections of Eq. Q on and yield [T] 

U^'^^e = - (e + P±) ^^U^' + (p^ - P\^ U.V^d^V^, (11) 

v^d^P\\ = - (P|| - p^_) d^v^ + {e + p±) v^f^df^ur (12) 

For one-dimensional motion along the beam axis it is sufficient to consider 
Eqs. (|10p . (|lip . and (|12p . Moreover, if the motion is boost-invariant, Eq. 
is satisfied automatically and we are left with Eqs. (fTUj) and ([TT]) only. 



3. Microscopic interpretation 

The anisotropic structure of the energy-momentum tensor appears nat- 
urally if partons are described by the phase space distribution function 



/ = /K^,^ , (13) 




where Aj_ and Ay are two different space-time dependent scales that may be 
interpreted as transverse and longitudinal temperature, respectively. The 
form p3p is valid in the local rest frame of the fluid element. The manifestly 
covariant form of / is 

/ = (14) 



^ The case E = was analyzed in Refs. [141 115] . 
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We note that as a special case of ()14p we obtain the formahsm of transverse 
hydrodynamics [IllITT], see Sect. [12] for more details. 

If the distribution function is known, one can calculate easily the particle 
density 



(27r)3-' \^A^' A| 
Similarly, we calculate the energy density 



JP± \P\\\\ (-,f,^ 



where Ep = ^V±_ + P\\ i ^-iid the transverse and longitudinal pressures 

"-fi'^J^], (17) 



(27r)3 2Ep ' I Ax ' All 



(27r)3 Ep' I Ax' A|| 



(18) 



Equations (fT5]) -p ^ indicate that we may use A_l and Ay alternatively 
with P± and Py. In Ref. [15] we have shown, however, that the two most 
convenient thermodynamics-like parameters are the particle density n and 
the anisotropy parameter x defined by the equation □ 



(19) 



With the help of those two variables we may express the energy density and 
pressures in a compact form |14^ [T5] , namely 




e(x, n) 
P|| (x, n) 
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R{x) 



+ xR'{x) 
- 2xR'{x) 



(20) 
(21) 
(22) 



Since the entropy density is proportional to the particle number density, an alternative 
choice is to use the entropy density a and x, see Eqs. (|3i p - (|33|l . 
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where the function R{x) is defined by the integral 

R{x) = x-y'J (23) 
and 3 is a constant defined by the expression 

Here = p±/X_l and = P\\/\\ are two dimensionless parameters. The 
range of the integration over ^_|_ and .^y is always between and infinity. 
The function R'{x) in Eqs. (j2ip and (|22p is the derivative of R{x) with 
respect to x. It is expected that R'{x) vanishes at x = 1, in this case the 
value X = 1 corresponds to equilibrium where P± = P\\- 



4. Entropy source 



To solve Eqs. Q and (jlOp we have to assume a certain form of the 
entropy source, S = S(o", x). In Ref. [Tl we proposed the following ansatz 

S = So(cr, x) = — = ^ , (25) 

A±A|| Teq ^JX Tgq 

where Teq is a time-scale parameter. The expression on the right-hand-side 
of Eq. (f25|) has several appealing features: i) it is positive, as expected 
on the grounds of the second law of thermodynamics, ii) it has a correct 
dimension, iii) it vanishes in equilibrium, where x = 1, iv) it is symmetric 
with respect to the interchange of with Ay (consequently S does not 
change if x — )• 1/x), finally v) for small deviations from equilibrium, where 
|x — 1| <C 1, we find 

fx - 11^ 

So(^) ^ ^7-^^- (26) 

'eq 

The quadratic dependence displayed in (j26p is characteristic for the 2nd 
order viscous hydrodynamics |j. Below we consider also the entropy source 
of the form 

{P±-P\\? a 
S = E(a,x)= ^ p p"^ ^, (27) 

where feq is another time-scale parameter that may be different from T^q. 

In addition to expressions ([25]) and (f27|l we also identify the entropy 
source in the model that has been recently proposed by Martinez and Strick- 
land [H]. This issue is discussed in Sect. 16.31 



* For the 2nd order viscous hydrodynamics, the entropy production is proportional to 
the viscous stress squared, E oc H'^ [T5]. On the other hand, for small x — 1 one can 
find that 11 is proportional to a; — 1 [12] , See also Ref. [TS] 
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5. Standard formulation and corrections for quantum statistics 



5.1. Anisotropic Boltzmann distribution 

Before we switch to the discussion of the effects connected with the use 
of different forms of the entropy source, it is interesting to discuss differ- 
ent forms of the distribution function which motivates the structure of the 
energy-momentum tensor used in our approach. Our previous calculations 
have been based on the following exponential form of the distribution func- 
tion 

f = go^w{-^Jvi/W+p\l^\)■ (28) 
Here qq is the degeneracy factor connected with internal quantum numbers 
such as spin, flavor, or color. Since we assume that the system formed at 
the early stages of collisions consists mainly of gluons, we take go = 16. 
Consistently with (j28p . the entropy density a has been obtained from the 
Boltzmann non-equilibrium definition 



gpAj A|| r f{i^A\) 



1 - In 



90 



(29) 



Comparing ()29p with the particle density obtained from (jlSp one gets 

o- = 4n. (30) 

In this case the parameter g obtained from Eq. (|24p equals go/ir'^ and 
Eqs. (|20p - ()22p may be rewritten in the equivalent form as 



4/3 



/ 2 \4/3 



P||(x,ct) 



\^9o I 



^ + xR'oi^) 
Ro{x) 



2xRq{x) 



(31) 
(32) 
(33) 



where Ro{x) is the function obtained from Eq. (|23p with the distribution 
(1281) . namely 



Roix) 



3gox 3 
27r2 



1 + 



X arctan ^/x — 1 
Vx^ 



(34) 



The function Ro{x) and its derivative are shown in Fig. [H we note that 
Ro{l) = Sgo/n'^ and, as expected, Rq{1) = 0. 
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Fig. 1. Function Ro{x) defined by Eq. (p4)) (upper part) and its derivative R'oix) 
(lower part). 



5.2. Anisotropic Bose-Einstein distribution 

The function (|28p may be considered as a Boltzmann equilibrium distri- 
bution that has been stretched (or squeezed) in the longitudinal direction. 
It is interesting to consider, in the similar way, the stretched (or squeezed) 
Bose-Einstein distribution. 



r 90 

I BE = — 



exp 



v/piAi+pJ/Aj 



The corresponding definition of the entropy density for bosons is 
5oAi A 



a 



27r2 



go go \ go J \ go 



(35) 



(36) 
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In this case the connection between the entropy density and the particle 
density is the same as in the isotropic gas of massless bosons, 



27r^ 



n. 



45C(3) 

Substituting Eq. ([35]) into Eqs. and ^ we find, see Sect. HSl 



gBE = C(3)-^, 



RBEix) = m Roix) 

Thus, for quantum statistics we obtain 



90 



Roix). 



e{x, a) 



Pi (x,a) 



B\{x,a) 



'45 C^/^(4) 



'45 C^/^{4)a 



4/3 

I Roix 

4/3 , 



'45 C^/4(4) 



4/3 



Roix) 



Roix) 



+ xR'Qix) 



2xRQix) 



(37) 

(38) 
(39) 

(40) 
(41) 
(42) 



We note that the numerical factor 45 C^''^(4)/(27r^) ~ 2.42 appearing in 
Eqs. (fl0|) ~ ([l2]l is very close to the factor 7r^/4 2.47 appearing in Eqs. 
([3T|l - (f33|) . hence, the equations defining the energy density and pressures in 
terms of the entropy density and x for the classical and quantum systems 
are practically the same. 



6. One-dimensional boost-invariant motion 

In the case of the longitudinal and boost-invariant evolution, the four- 
vectors and may be defined by the expressions 



= (cosh T], 0, 0, sinh rj), 
= (sinh ?7, 0, 0, cosh r/). 



where rj is the space-time rapidity. 



1, t + z 

r? = — in . 

' 2 t-z 



(43) 
(44) 



(45) 
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Substituting Eqs. (gS]) and (03]) in Eqs. ([9]) and 1^ we obtain 

— = ^, 46 

or r 

da a ^ 

— + - = S, 47 
dr T 

where r is the proper time 



r = (48) 

6'. i. Standard formulation 
Substituting Eqs. ([31]) and ([Ml) in (06]), and using Eq. ([25]) for the 



entropy source, we obtain one ordinary differential equation for x, 

dx 2x 4:Ho{x) 



dr T 3rcq 



(49) 



We have introduced here the function Hq{x) defined by the expression 

i?o(x)S Rojx) (1 - 

For small deviations from equilibrium we may use the expansion 

Ho{x)^-{x-l) + —{x-lf + .... (51) 
Equation (|49p was analyzed in detail in Ref. [1], where we showed that 

X ^ 1 for r ^> Teq. 

6.2. Modified entropy source 

Let us now turn to the discussion of the entropy source defined by the 
difference of pressures. If we use Eq. ([27D instead of Eq. ([2^ . we obtain 

dx 2x 4:H(x) . , 

(52) 



dr T 3feq 



where 



fr(. ^ 9x'Ro{x)R',{x) 

{Ro{x)/3 + xR'o{x)){Ro{x)/3-2xR'o{x))- ^ ' 
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-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

^ = x-l 



Fig. 2. Function Ho defined by Eq. ([5(I| (solid line) compared to the rescaled 
functions H and Hms defined by Eqs. (j53p and (1501) (dashed and dotted lines, 
respectively). The arguments of the functions have been shifted to ^ = x — 1. 



Similarly to ([51]) we may expand (|53|1 around x = 1, 

H[x)^f{x-l) + ^{x-lf + ---. (54) 

The requirement that the two definitions of the entropy source lead to the 
same behavior in the region close to equilibrium leads to the condition 



64 
25 



Teq. (55) 



6.3. Martinez and Strickland model 

A similar model to ADHYDRO has been introduced recently by Mar- 
tinez and Strickland (MS) [12^ JjiJ. This model is based on the consideration 
of the two first moment of the Boltzmann kinetic equation with the colli- 
sion term treated in the relaxation time approximation. The inverse of the 
relaxation time is denoted by F. 

In this approach the evolution of matter is described by two coupled or- 
dinary differential equations for the momentum anisotropy ^ and the typical 
hard momentum scale Phard- They may be written in the compact form as 

d^ = —r 3 ' ^''^ 

1 ^^'hard_ri?;,3(0^^^^^^_ ^^^^ 



Phard 
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Fig. 3. Different entropy source terms normalized by (j/tcc^. The relaxation times 
are connected by the formula ([5^ . 



The parameter ^ is related to x, 

i = x-l. (58) 
Similarly, the function i?Ms(0 is simply connected to Ro{x), 



RusiO = \ 



1 arctan JT 
+ 



1 + C VI 



;^(l + 0-'/'^o(l + e). (59) 



In order to show similarities with our approach, in Eqs. ()56p and ()57p we 
have introduced the function Hyis that is an analog of Hq{x), 

3(e+i)(ve+T<s(0-i)^Ms(e) ^ ^ 
"""^^^^ = — 2i?Ms(e) + 3(e + i)i?;.s(0 — • ^''^ 

Close to equilibrium, the function Hms{£.{x)) has the following expansion 
HusiCix)) ^l{x-l) + ^{x-lf + ---. (61) 



Thus, comparing Eqs. (jMj) and (j56|l we conclude that our standard ap- 
proach is equivalent to MS (near equilibrium) if 

" (62) 



2Teq 5Teq 
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The plots of the functions Hq{x{S^)), H{x{$,)) and Hms{x{(,)), rescaled by 
the factors defined in (j62p . are shown in Fig. [2j As it has been aheady 
discussed in [2] the essential difference between these functions is observed 
only for very large anisotropics. 

It is instructive to identify the entropy source in the Martinez and Strick- 
land model. The entropy density in this model has the form a = ^Phard^^^^^- 
With the help of Eqs. (|56|) and (|57|) we reproduce Eq. (|47|) with 



^MS 



RusiO 3(1 + 



Thus, comparing Eqs. ([SH]) and ([^5]) we find that the Martinez and Strick- 
land model fits exactly to the ADHYDRO scheme if the entropy source is 
defined by ([63|1 . 

Moreover, expanding functions RusiO ^^'^ Hms{S,) for small ^, and 
using Eq. ([62]) . we find that Ems agrees with Sq and S up to quadratic 
terms in ^. Different entropy source terms normalized by cr/req are shown 
in Fig. [3l Once again, the relaxation times are identified with the help of 
Eq. 



7. Results for the boost-invariant case 

In this Section we present numerical solutions of Eqs. (j49p . (j52p and (j56p . 
As discussed above, these equations correspond to three different forms of 
the entropy source defined by Eqs. ([2^ . ([TTp . and respectively. The 
parameter Teq is set equal to 0.25 fm. Using this value we find that the 
system equilibrates at the proper time of about 1 fm. The values of the 
time-scale parameter fgq and the inverse relaxation time F are connected 
with Teq by Eq. ([62]) . 

In all the cases the initial conditions are specified at the initial proper 
time To = 0.2 fm. We assume that the transverse pressure dominates over 
the longitudinal pressure during the very early stages of the evolution of 
matter and take x{to) = xq = 100. In the upper part of Fig. [J] we show the 
time dependence of x. The solid line describes the standard ADHYDRO 
result, with S = Sq, the dashed line describes the result obtained with the 
modified entropy source in ADHYDRO, with S = S, and the dotted line 
describes the time evolution of x in the MS model where S = Sms- 

The middle part of Fig. S] shows the time evolution of the ratio of the 
longitudinal and transverse pressure, and the lower part of Fig. U] shows 
the time evolution of the entropy density normalized to the initial entropy 
density. The notation is the same as in the upper part of Fig. [H 
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Fig. 4. Time dependence of the asymmetry parameter x (upper part), ratio of the 
longitudinal and transverse pressures (middle part), and entropy density normal- 
ized to the values obtained from the Bjorken model, tiBj = ctoTq/t (lower part). 
The standard ADHYDRO calculations are shown as solid lines, the results of AD- 
HYDRO with the modified entropy source are represented by dashed lines, and the 
results for the MS model are shown as dotted lines. 
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We observe that the time dependence of x is very much similar in the 
three cases, suggesting that the particular form of the entropy source has 
small effect on the time evolution (note the double logarithmic scale). This 
behavior is caused by the initial large entropy production that leads to a 
fast decrease of x. Later, when x becomes closer to unity, all the schemes 
are equivalent and imply a very similar time development. 

We find, however, that the equilibration of pressures is a little bit delayed 
in the MS model, as compared to the standard version, and this is reflected 
in a larger entropy production (relative to the initial entropy density). On 
the other hand, the equilibration in the modified model where S = E is a 
little bit accelerated. 



8. Non-boost-invariant evolution 

In the case of purely one-dimensional, non-boost-invariant motion we 
define 

U^' = (cosht?(r,??),0,0,sinh??(r,r?)) , (64) 
= (sinhi?(T, r/), 0, 0, cosh i?(r, r/)) , (65) 

where ??(r, rj) is the fluid rapidity which depends on proper time and space- 
time rapidity. Equations ([5^ and (|65|) guarantee that Eqs. (HD are auto- 
matically satisfied. 

Using Eqs. and ([U^ in Eqs. ^ and (fTUj) we derive three partial 
differential equations for (T(r, rj), 7?(r, 77), and a;(r, if). The structure of those 
equations has been discussed recently in Ref. ^ and in this paper we skip 
details describing their derivation and form. It is important, however, to 
present here our initial conditions. 

The initial conditions for a non-boost-invariant evolution are defined by 
the three functions of space-time rapidity: o"(ro,r/), ?9(to,??), and x(ro,r/). 
The initial entropy profile is taken in the form [111 [2T] 



(j(to,?7) = fToexp 



(66) 



where 9 is the step function, the parameter A?] defines the half width of 
the initial plateau in space-time rapidity, and drj defines the half width of 
the Gaussian tails on both sides of the plateau. To match the rapidity 
distribution measured by BRAHMS [22j we use the values: At/ = 1 and 
5r] = 1.3. The value of the initial central entropy density is obtained from 
the condition that the initial energy density is 100 GeV/fm^, 

/ 2 \4/3 

eo = 100GeV/fm3 = ^ R{xo). (67) 

V / 
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Fig. 5. Upper parts; time evolution of the anisotropy parameter x for two different 
values of the space-time rapidity rj = (left) and 77 = 6 (right). Lower parts: time 
dependence of the ratio of the longitudinal and transverse pressures, P\\/P±, again 
for 77 = (left) and 77 = 6 (right). The solid lines represent the standard calculation 
where E = Eq, the dashed lines show the results with the modified entropy source 
where E = E, and the dotted lines show the results for E = Ems- 



Here, xq is the initial value of the anisotropy parameter at 77 = 0. In this 
Section we assume that the initial profile of x is constant, 

x(ro,r/) = xo = 100. (68) 

The time-scale parameter Tgq has the same value as in the previous Section, 
and the initial fluid rapidity profile is taken in the form 

^{to,v)=V- (69) 

The results of our numerical calculations done with three different forms 
of the entropy source are shown in Figs. [5]and[6j The solid lines represent 
the standard calculation where S = Sq, the dashed lines show the results 
with the modified entropy source where S = S, and the dotted lines show 
the results for the case S = Sms • In the two upper parts of Fig. [S] we 
show the time dependence of the anisotropy parameter x for two different 
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Fig. 6. Time evolution of the entropy (upper part) and fluid rapidity (lower part) 
profiles in rj. Similarly as in Fig. O the solid, dashed, and dotted lines correspond 
to three different forms of the entropy source used in the numerical calculations. 

values of the space-time rapidity, rj = (left) and rj = 6 (right). We observe 
that this dependence is very much similar to the dependence obtained in 
the boost-invariant calculation, see the upper part of Fig. HI Thus the 
rapidity dependence of the initial entropy profile does not affect the time 
dependence of x in a noticeable way. In the two lower parts of Fig. [5] 
we show the time dependence of the ratio of pressures, P\\/Pj_- Again, we 
observe weak rapidity dependence, compare the middle part of Fig. [H 

In Fig. [6] we show the entropy density a and the fluid rapidity ■!? as 
functions of space-time rapidity for the initial proper time r = tq, a bit 
later proper time r = 1.2 tq, and r = 3ro (for o") or r = 10 fm (for ■& — rj). 
At T = 1.2 Tq the entropy profiles are higher than the initial profiles — at 
such an early time the entropy is produced in the system. For r = 3 tq the 
entropy profiles lie below the initial profiles — the entropy is not produced in 
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Fig. 7. Time evolution of the anisotropy parameter x (upper part) and the ratio of 
pressures (lower part) for three different values of the space-time rapidity: = 
(solid lines), rj = 4.9 (dashed lines), and -q = & (dotted lines). 

the system anymore, and the entropy density drops due to the longitudinal 
expansion of the system. We observe that the results obtained with the 
three forms of S are very much similar, indicating that the specific form 
of the entropy form has small effect on the evolution of matter also for the 
non-boost-invariant systems. 

9. Rapidity dependent anisotropy 

In the previous Section we have analyzed the case where the initial 
anisotropy parameter is independent of space-time rapidity, x{To,r]) = xq. 
In this Section we consider the case where the initial anisotropy profile 
depends on rj but other initial conditions are exactly the same. We assume 





Teq=0.25fm 




;? = 




rj = 4.9 




'/ = 6 _ 
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that the initial rj profile of x is the same as the initial entropy density profile, 
namely 

x(To,r]) = xo , (70) 

where cj(ro,r/) is given by Eq. (j66|) . With xq = 100 Eq. ([70]) implies that 
the initial anisotropy is large at the center of the system (central region) 
and small at the edges of the system (fragmentation regions) . 

In Fig. [7|we show the time evolution of the anisotropy parameter x and 
the ratio of pressures for three different values of the space-time rapidity: 
rj = (solid lines), rj = 4.9 (dashed lines), and rj = 6 (dotted lines). The 
behavior of the system in the central region (rj = 0) is very much similar to 
the behavior of the boost-invariant systems studied before. On the other 
hand, at large rapidities the time evolution of the physical quantities is quite 
different, which is caused by a different starting value of x. For small initial 
values of x we observe first the growth of x with time (above unity) and 
only later a decrease towards unity. Anyway, we observe that at r ~ 0.5 fm 
the values of x become practically independent of r]. 

In Fig. [8] we show the time evolution of the entropy and rapidity profiles 
for the initial conditions specified by Eqs. ([66]) . ([68l) . and ([69|) (solid lines) 
and by Eqs. ([66]) . ([69]) and ([70]) (dashed lines). In the two considered cases 
the entropy profiles behave in the very similar way. The behavior of the 
fluid rapidity is different only at the beginning of the evolution for large 
rapidities. 



10. Conclusions 

In this paper we have discussed the assumptions of the recently for- 
mulated model that may be used to describe the evolution of a highly- 
anisotropic fluid formed at the very early stages of relativistic heavy-ion 
collisions. We have considered three forms of the entropy source which de- 
fines the dynamics of the system. We have found that the numerical results 
depend very weakly on the particular form of the entropy source. 

We have also studied microscopic interpretation of the model, discussing 
different forms of the microscopic phase-space distribution function. We 
have found that corrections for quantum effects lead to minor quantitative 
differences as compared to the standard case where the phase space distri- 
bution is described by the (stretched or squeezed) Boltzmann distribution. 

We have performed the numerical calculations for the cases where the 
initial anisotropy depends on space-time rapidity. In such cases the evolu- 
tion of the system is more complicated but the tendency of the system to 
equilibrate is not changed. 
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Fig. 8. Time evolution of the entropy and rapidity profiles for the initial conditions 
specified by Eqs. ([SS]), (|M1), and dSS]) (sohd hnes) and by Eqs. and ([70)1 

(dashed lines). 

Altogether, we have found that the results obtained with the ADHY- 
DRO model are stable against various modifications of the model compo- 
nents such as: the entropy source, the form of the microscopic phase-space 
distribution function, or the form of the initial conditions. We consider such 
stability as an attractive physical feature of the model proposed in [T] . 

11. Appendix 1 

In this Section we give the explicit form of the Lorentz transformation 
which leads from the center-of-mass frame of the colliding nuclei to the 
local-rest-frame of the fluid element. This transformation is constructed as 
a sequence of the longitudinal Lorentz boost, rotation around the z-axis, 
and the transverse Lorentz boost. In this way the complete transformation 



21 



does not change the longitudinal direction. 

The four- vectors C/'* and may be written in the following form 



U^" = cosh??iicosh7?i 1 



tanh ■!? I cos 6 tanh 'd i cos ( 



cosh i9 1 1 



coshi?ii 



■,tanh??|| 



V^" = cosh (tanh 0,0, l) 



(71) 



where ??|| is the longitudinal fluid rapidity and ■d± is the transverse fluid 
rapidity. 

The Lorentz transformation leading to the local rest frame consists of 
the Lorentz boost along the z-axis, 



7- M _ 



( coshi?|| 



\ — sinh-t?!! 



-sinht?|| 

1 
1 

cosht?,. 



the rotation around the 2;-axis, 



(72) 



{xy)v 



1 1 











\ 





COS (j) 


sin (j) 










— sin^ 


COS (p 







V 








1 





(73) 



and the Lorentz boost along the x-axis 



( cosh 
— sinh i?j 


V 



-sinht?_L \ 
cosht?_L 








By direct multiplication one may check that 



1 
1 ) 



L 



(x)l3 ^(xy)-i {z)ii 
)/3 



(1,0,0,0), 



(74) 



(75) 



In order to move from the local rest frame to the center-of-mass frame we 
perform the inverse transormation. 



12. Appendix 2 

The concept of ADHYDRO includes, as a special case, the framework 
of transverse hydrodynamics. In the latter case, we assume that the distri- 



22 



bution function ()14p has the following factorized form 



(76) 

Here uq is a normalization constant and S is the Dirac delta function. Equa- 
tion ()76p allows us to make the identification 

All = 1. (77) 
Using our general definition of the function R{x), see Eq. ([23|) . we find 

R{x) = rxl/^ xi?'(x) = (78) 

6 



where r is a constant. Equations (|78|1 used in (|20|) ~ (|22|) immediately yield 

^± = |, ^'11=0. (79) 

Hence, the ansatz (j76p corresponds to the situation where the longitudinal 
pressure vanishes and the system may be treated as a collection of indepen- 
dent clusters which expand transversally. 

The parameter A_l describes temperature of two-dimensional transverse 
clusters. For the purely longitudinal expansion, the parton density drops 
like l/r. Since n ~ -^i-^H = ^'±/'^^ the value of A_l remains unchanged. 
Clearly, only transverse expansion can lower cluster's temperature. More- 
over, in such a boost-invariant case x ~ r^, hence the ratio of the transverse 
and longitudinal pressure always grows and the system cannot reach stable 
isotropic equilibrium. 



13. Appendix 3 

In this Appendix we discuss consequences of replacing the classical distri- 
bution function by its quantum analogs, the Bose-Einstein and Fermi-Dirac 
distributions. The integrals over the distribution function (|35p may be done 
by expanding ([35]) in a geometric series, 

oo 

fBE = 50 E (-ny^pi/Ai+p|/Aj) , (80) 

n=l 
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and by integrating ()80p term by term. For example, starting from the 
general definition of the function R{x), see Eq. ([23|) . we obtain 



n=l 

(81) 

Introducing the new integration variables: x± = n^j_ and X|| = n^y we find 

n=l 

(82) 



Rbe{x) = gox ^^^^^ J ^^^^-^^^W y^a;| + xx^ exp [^-^x^ + x| 



The sum over n factorizes and gives the Riemann zeta function C(4)- As a 
consequence, we obtain the simple relation to the "classical" result, 

Rbe{x) = am^ix) = '^Mx). (83) 

We note that a similar technique may be used in the calculations involv- 
ing the Fermi-Dirac distribution. In this case we have 

JFD 



n=l 

and Eq. ([23]) leads us to the expression 



exp(y^pi/Ai+p|/Aj)+l 

oo 

50 E(-l)"^'e^P (-n^pi/Ai+^^f/Aj) , (84) 



RMx) = mx-'/'f^i-iT'-' 

n=l 

(85) 

As in the standard treatment of the Fermi-Dirac distributions in statistical 
mechanics, we may change the summation in the following way [20j 

E (-ir^--- = E E ••• 

n=l,2,3,... n=l,3,5,... n=2,4,6,... 

= E •••-2 E •••• 

n=l,2,3,... n=2,4,6,... 



24 



In the last sum we change n to 2n and obtain 

n=l 



n=l' 



(87) 



Now, in the first integral on the right-hand-side of (|87|) we change the inte- 
gration variables to x± = n^_L and = i^i t^is second integral we 
change to x± = 2nS^± and xy = 2n^||. In this way we find 



Rfd{x) 



n=l 



27r2 



where the factor of 7/8 follows from the term 1 — 2^"^, where p is the total 
power of ^j^s and .^y's appearing in the second integral in ([87|l . Thus, the 
final result is 

7C(4) „ _ 7^4 



Rfd{x) 



-Ro{x) 



720 



i?o(x). 



(89) 



We stress that the definition of Ro{x) involves the number of the internal 
degrees of freedom hence when dealing with both bosons and fermions 
at the same time, one should introduce Ro{x) separately for bosons and 
fermions and add them together. 



Rbe+fd{x) 



vr 
60 



1 + 



: arctan ^/a 



1 



1 



(90) 



Here and g^^ are the numbers of internal degrees of freedom of bosons 
and fermions, respectively. 

The parameter g defined by Eq. ()24p may be also calculated for the 
case of massless bosons and fermions. Following the same steps as in the 
calculation of the functions Rbe{x) and Rfd{x) we obtain 



gBE 

QBE = C(3) — 2~' 



9FD 



C(3) (] 



-2\ % 



(91) 



(92) 
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and for the sum 



9BE+FD = 



C(3) 



( 




,FD 



) 



(93) 



Equations ([90]) and ([93]) may be used for R{x) and g appearing in Eqs. (f20l) - 
()22|) if the system of quarks and gluons with the same anisotropy parameter 
X and the total density n is described. 
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